Asymptotic analysis of mode-coupling theory of active nonlinear microrheology 
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We discuss a schematic model of mode-coupling theory for force-driven active nonlinear microrhe- 
ology, where a single probe particle is pulled by a constant external force through a dense host 
medium. The model exhibits both a glass transition for the host, and a force-induced delocal- 
ization transition, where an initially localized probe inside the glassy host attains a nonvanishing 
steady-state velocity by locally melting the glass. Asymptotic expressions for the transient density 
correlation functions of the schematic model are derived, valid close to the transition points. There 
appear several nontrivial time scales relevant for the decay laws of the correlators. For the nonlin- 
ear friction coefficient of the probe, the asymptotic expressions cause various regimes of power-law 
variation with the external force, and two-parameter scaling laws. 

PACS numbers: 82.70.-y 64.70.pv 83.10.-y 



I. INTRODUCTION 

Microrheology is a modern technique that allows to 
probe complex fluids on mesoscopic length scales. One 
inserts a probe particle, typically /im-sized, into a host 
liquid of constituents that are roughly of the same size 
(such as colloidal dispersions or biophysical fluids) . Mon- 
itoring the motion of the probe, one can infer local visco- 
elastic response functions of the host liquid. A partic- 
ularly compelling extension of the technique is called 
active microrheology: here, the probe is subjected to a 
controlled external drive. This is most conveniently ap- 
plied by using laser tweezers or magnetically susceptible 
probe particles, or even by tailoring colloidal probes such 
that they undergo self-driven motion due to chemical pro- 
cesses [IH3]. Active microrheology has become a major 
tool in biophysics and for colloidal model systems [IH^- 

Here we focus on force-driven active microrheology, 
where a constant external force F_^y. is applied to the 
probe particle. A natural quantity to observe then is the 
resulting probe velocity v_, and specifically its ensemble- 
averaged stationary value, {v)t^oo- The influence of the 
host liquid is characterized by a friction coefficient, 

av)t^oo = Zex ■ (1) 

For typical soft-matter systems, thermal fluctuations give 
rise to forces in the range of pN; it is hence easy to drive 
the system into the nonlinear nonequilibrium regime us- 
ing active microrheology. In Eq. ([T]), the friction coeffi- 
cient C(-Fex) then becomes a function of the applied force. 
This makes analysis of the experiment vastly more diffi- 
cult, since one needs to employ a theory of the nonlinear 
probe-bath interactions. If that is available, the tech- 
nique can, however, be rewarding, as it gives access to 
much more information about the complex host liquid 
than a linear-response setup would. Also, compared to 
other nonlinear-response techniques, such as macroscopic 



rheology, one gains access to the microscopic mechanisms 
relevant for the dynamics [8.. 

Recently, a microscopic theory for the force-driven ac- 
tive nonlinear microrheology has been proposed [9^, us- 
ing a combination of the integration-through transients 
(ITT) scheme together with approximations inspired by 
the mode-coupling theory for the glass transition (MCT) . 
ITT expresses the nonlinear friction coefficient C(-Pcx) via 
a relation of (generalized) Green-Kubo type, through a 
transient force autocorrelation function. This is a cor- 
relation function taken with the equilibrium ensemble, 
but the full nonequilibrium dynamics. Governed by the 
idea that in a dense system, structural relaxation via 
density fluctuations is the dominant slow dynamical pro- 
cess, MCT provides an approximate closure for this cor- 
relation function in terms of (again, transient) density 
correlators. These are calculated by a set of nonlinear 
integro-differential equations. As a result, the nontrivial 
relaxation pattern predicted for the density correlators 
directly gives rise to nonlincarities in the friction coeffi- 
cient. 

So far, the full MCT-ITT equations for active mi- 
crorheology have proven not amenable to extended nu- 
merical treatment. The strategy used in Ref. [9] there- 
fore was to implement an ad-hoc simplified model that 
reduces the complexity down to a single fluctuation mode 
in the hope of retaining all nontrivial mathematical fea- 
tures of the original set of equations. Introducing a few 
adjustable parameters to this so-called schematic model, 
a successful quantitative analysis of available computer- 
simulation and experimental data provides an ex poste- 
riori justification for doing so. The original schematic 
model proved too restrictive in certain aspects. Recently, 
an improved version of the model has been presented [TO] , 
which allows convincing fits of the data in all accessible 
regimes by taking into account some aspects of the force- 
induced spatial anisotropy in the dynamics. 

The data show a strong nonlinearity in ((Fcx)- in a 
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relatively narrow range of forces, close to the glass tran- 
sition the friction coefficient drops by orders of mag- 
nitude, separating a near-equilibrium low-force regime 
from a high-force regime. The schematic MCT models 
interpret the resulting strong drop as a precursor to a 
probe-delocalization transition: inside a glassy host, the 
probe particle is held in a nearest-neighbor cage (formally 
C — >■ oo at low forces, for the ideal-glass case) that can 
sustain a finite amount of external force. Once the ap- 
plied force exceeds a critical threshold F^^, the probe's 
nearest-neighbor cage is forced open so that a finite mean 
velocity results (and ( drops to a finite value). The 
threshold force is thus interpreted as a measure of the 
cage strength. Since cages are formed in a highly collec- 
tive process involving the host particles and their local 
structure, active nonlinear microrheology is in principle 
a unique tool to probe the local rigidity of that struc- 
ture. In fact, F^^ = 0{5QkBT/a) has been measured 9J, 
greatly exceeding the typical force scale of thermal fluc- 
tuations, k^T/a (here, cr is a typical host-particle size). 

In order to understand the asymptotic behavior of 
C(Fex) close to the critical force within ITT-MCT, the 
asymptotic behavior of the correlation functions close to 
this delocalization transition has to be understood. In 
essence, a two-parameter scaling prescription is sought 
for, since one is dealing with both the distance to the 
glass transition, and the distance to the delocalization 
transition, as small parameters. 

In this paper, we present an asymptotic treatment of 
the schematic ITT-MCT equations for force-driven ac- 
tive nonlinear rheology. We discuss the set of equations 
that has been solved numerically in Ref. [TU], where the 
model was shown to describe both computer simulation 
and experiment quantitatively. The calculation proceeds 
by considering the double limit of approaching the glass 
transition and the delocalization transition, which yields 
two small parameters and various asymptotic results de- 
pending on their ratio and sign. The techniques we use 
are similar to those that have been used earlier to de- 
rive different two-parameter scaling laws within MCT, 
e.g., for the extended mode-coupling theory including a 
schematic hopping term [TT], or for the ITT-MCT equa- 
tions describing the macrorheology for a given constant 
shear rate JT2]. 

Our analysis, however, has aspects that differ from the 
above-mentioned cases. Due to the nonequilibrium na- 
ture of the problem, the time-evolution operator is non- 
Hermitian, and gives rise to complex-valued correlation 
functions. Usually, in the regime of structural relaxation, 
one can safely assume the correlation functions appearing 
in the theory to be real-valued and completely monotone; 
this holds rigorously for overdamped short-time dynam- 
ics as applicable to colloidal suspensions in equilibrium 
[131 m]. These properties are used in the asymptotic 
expansion, for example to ensure that the singularities 
to be discussed belong to a certain class of bifurcations 
[15j . Even in the macro-rheology of colloidal suspensions, 
where external flow is represented by a non-Hermitian 



generalized Smoluchowski operator, taking into account 
the mechanism of shear advection separately allows one 
to return to real-valued correlation functions and to a 
scheme of asymptotic expansions that closely mirrors the 
one followed in equilibrium jT^18| . 

In the equations for active microrheology, some as- 
sumptions entering the standard discussion of MCT are 
no longer valid. Consequently, the mathematical classifi- 
cation of the transition between localized and delocalized 
probe particles (in the idealized glass) is still open. We 
restrict ourselves here to a certain schematic model that 
is inspired by, but not necessarily mathematically equiv- 
alent to, the microscopic equations presented in Ref. [S]- 
Since the model has been successfully used in data anal- 
ysis, the restriction appears plausible. 

The paper is organized as follows: in Sec. |llj we sum- 
marize the equations defining the model. The long- 
time limits of the correlation functions, characterizing 



the glassy and localized states, are discussed in Sec. Ill 



Sections pV| and |V| are devoted to deriving asymptotic ex- 
pressions for times large compared to the single-particle 
relaxation time, valid on intermediate- respectively long- 
time windows that open upon approaching the transition 
points. These are the analogues to the common MCT 
scaling laws referred to as [3- and a-scaling. In Sec. |VI| 
we transfer these results to a two-parameter scaling law 
for the friction coefficient, after which Sec . | VII| concludes . 



II. SCHEMATIC MODE-COUPLING THEORY 

We summarize the main equations defining the 
schematic MCT model for active nonlinear microrheol- 
ogy. Following a generic integration-through-transients 
(ITT) scheme and the notion that the slow dynamics in 
the vicinity of the (colloidal) glass transition is domi- 
nated by density fluctuations [TS], the central quantities 
of the model are the transient density correlation func- 
tions. While in the liquid state, these eventually decay 
to zero, in the glass they attain a flnite positive long-time 
limit called the glass form factor or nonergodicity param- 
eter /. Even in the liquid, the correlation functions stay 
close to its value fc at the transition over an increasingly 
large time window as one approaches the transition. The 
quantity / and its tagged-particle counterpart will 
play a central role in the discussion to follow. 

In the speciflc model we choose, there is one correla- 
tor (j){t) mimicking the dynamics of the host liquid and 
determining the glass form factor / = limt_>.oo 0(i). Fur- 
thermore, this dynamics is taken as identical to the equi- 
librium one. Clearly, in the thermodynamic limit and 
assuming that the probe-host interactions remain suffi- 
ciently short-ranged, the host-liquid dynamics is in the 
ensemble average unperturbed by the external force that 
is applied to the probe particle only. The macroscopic 
state of the host (liquid or glassy) will then be deter- 
mined by the equilibrium coupling coefficients only. The 
equation of motion for the host-liquid correlator takes 



3 



the form 

dt(l>{t) + r |</)(i) + ^ m{t - t')dt'(l)(t') dt'^ = . (2a) 

Equations of this type can be derived from microscopic 
starting points using a Mori-Zwanzig projection operator 
scheme [TS], where the slow relaxation is modeled by a 
memory kernel m{t) determined by the fluctuating forces. 
The coefRcient F is a relaxation rate of the short-time 
dynamics, and will be set to unity in the calculations 
below (thus defining the unit of time). The central idea of 
MCT is to approximate the memory kernel by a bilinear 
form of the correlators themselves, expressing the notion 
that a slow decay of density fluctuations leads to and 
hinges upon slow decorrelation of fluctuating forces. In 
our discussion we adopt the so-called F12 model, 

m{t) ^ vi(}){t) + V2(l){tf , (2b) 

with positive parameters (wi,W2). With this choice, the 
full range of asymptotic behavior expected from the mi- 
croscopic MCT close to ordinary liquid-glass transitions 
is reproduced [15] . The use of one single mode to describe 
the host dynamics embodies our assumption that the sys- 
tem remains homogeneous and isotropic in the ensemble 
average. 

Spatial isotropy clearly breaks down for the motion of 
the probe particle. Therefore, in the schematic model of 
Ref. [10], two tagged-particle correlation functions <j)%{t) 
were introduced, labeled by a G {|j , _L} to indicate the 
separate role played by fluctuations in direction of, and 
perpendicular to, the applied force. The equations of 
motion then read 

9tC(t) + Lj^ |c(i) + J^mUt- t')dt4l.{t') rfi'l = . 

(3a) 

Here, the initial relaxation rates are (jJ±_ — V g and wy = 
rs(l — iK||Fox), which fixes the unit of energy as k-gT = 1. 
We furthermore set Fs = 1 for simplicity, suitable for a 
probe that is almost identical to the host particles. K|| 
adjusts the scale of forces entering the schematic model. 
The slow dynamics of the probe arises as a consequence 
of slow dynamics in the host liquid, and is modeled by 
memory kernels containing a linear coupling to (/'(i), 

m|(t) - (vl<j,f{t) + v'2^l{t)) 0(t)/(l - iK\\F,^) , (3b) 

ml{t) = {vtcf>l{t) + vIRc^it)) (/.(i)/(l + (k^F,^)^) . 

(3c) 

The specific choice of terms entering Eqs. (|3| is rooted 
in symmetry considerations based on the full microscopic 
MCT model Specifically, for F^^ = 0, the model 

reduces to a well-studied schematic model of tagged- 
particle motion close to the glass transition, the Sjogren 
model [inj. In this case (^y(t) = <^j_(i), both are real- 
valued, and only a single coupling strength v'^ = vf + w| 



is relevant. Assuming that the probe motion retains rota- 
tional symmetry around the axis set by -Fox in the ensem- 
ble average, it is seen that (f'±{i) remains a real- valued 
function. On the other hand, a nonvanishing net dis- 
placement of the probe along the direction of the force 
is expected, which results in the modulation of the cor- 
responding Fourier-transformed density-fluctuation cor- 
relations with a complex phase. Hence, (t) will be 
complex-valued. Furthermore, Eqs. ^ obey the ex- 
pected symmetries under reversal of the applied force, 
Fox I— ^ —Fcx- The parameters Ky and k± characterize the 
forces relevant for inducing the decay of fluctuations in 
the two directions. They are introduced to achieve quan- 
titative fits of simulation and experimental data [10] . We 
fix Kj_ = 1/2 and K|| = 1 for the numerical calculations 
presented below. 

We are thus left with a three-correlator model where 
four parameters enter that are relevant for an under- 
standing of the qualitative long-time behavior: two of 
them, {vi,V2) model the approach to the glass transition 
in the host liquid. In the parameter space of the F12 
model, a line of ideal glass transitions (wi,c,i'2,c) exists, 
and it is useful to characterize the distance to any cho- 
sen point on this line by a single distance parameter e. 
For the discussion below, let us fix V2,c = 2 implying 
vi,c = 2{V2-1) [H], and set {vi,V2) = (wi,c, W2,c)(l 
As usual, e < signals liquid states, while e > holds 
for glassy states. We will denote quantities calculated at 
glass-transition points of the host-liquid model by sub- 
scripts c. 

The other two parameters, {vl,V2), represent the cou- 
pling of the probe to the host liquid; they will, among 
other things, also reflect a non-trivial size ratio between 
probe and host particles. We employ the simplification 
that was used in Ref. 10 and set = 2 for numerical 

calculations. Note that setting u| = in Eqs. ^ reduces 
the model to the one originally proposed in Ref. [9], not 
taking into account the role of probe-density fiuctuations 
perpendicular to the force direction. 

While the transient correlation functions can in prin- 
ciple be measured and have been evaluated in computer 
simulation [SI [TO] , an experimentally more easily acces- 
sible quantity is the friction coefficient (. A straightfor- 
ward adaption of the microscopic expression, used in }lO] , 
is C = 1 + AC with 

Jo 

poo 

+ (l-/i) / (l){t)Re(Pl{t)dt. (4) 
Jo 

This uses that in our choice of units, the solvent friction 
experienced by the free particle is unity, C = 1- In the mi- 
croscopic theory, is given as an angular average over a 
force-autocorrelation function; approximating the latter 
in terms of density-pair modes, one gets a Fourier-space 
integral over anisotropic coupling coefficients. In writing 
Eq. Q, we assume that this integral is qualitatively dom- 
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FIG. 1. Schematic-model correlation functions <j)'^^{t) and 
<l>±{t) for host-liquid parameters close to and at the glass 
transition, e = — 10~^ (solid lines) and e = (dashed), for 
probe-coupling coefficients v" = 30, k,± — 0.5, and «:|| — 1. 
For (^jj(t), both real and imaginary part are shown in sep- 
arate panels. Curves in the order of decreasing relaxation 
time correspond to Fa = 0, 1, F^^ — 0.1, F^^, 12, and 80, 
where = 6.5735. In the panels showing the real parts, the 
host-liquid correlation functions 0(t) are added for the liquid 
(dash-dotted lines) and glassy (dotted) state. 



inated by contributions from the two modes considered 
in the schematic model, and that in particular the quali- 
tative features of the correlators a —\\ and a —J- are not 
restricted to zero-measure portions of wave- vector space. 
The parameter fi allows to reweight these contributions, 
which does not qualitatively change the features close to 
the delocalization transition. Following Ref. [10] we set 
/i = 1/2 in our calculations. 

Figure [l] displays exemplary correlation functions of 
the schematic model both in the liquid state (solid lines) 
and at the liquid-glass transition point (dashed). Nu- 
merical solutions of Eqs. ([2]) to ([3| are obtained by inte- 
grating in the time domain, using a repeated doubling of 
the integration step to allow covering a large number of 
decades in time. The algorithm is a straightforward gen- 
eralization of the one used in previous MCT calculations 



Considering first F^x = 0, the correlation functions 
are all real, and Fig. [T] demonstrates the two-step decay 
typical for glassy structural relaxation: at times much 
larger than those associated with single-particle motion, 
t ^ 1/r, a window of structural relaxation opens. Corre- 
lators first decay towards a finite plateau, identifying by 
4>%{i ~ ia) ~ fa the so-called /3-relaxation regime. The 
time scale t„ diverges approaching the glass transition. 
In the glass, and right at the glass transition (e = 0), the 
correlation functions never decay from their plateau. In 
the liquid, a final decay to zero sets in on time scales large 
compared to those of the (5 relaxation, t ^ t^. This iden- 
tifies the Qf-relaxation window t/t'^ = 0(1) and a second 
time scale t'^ that diverges faster than to- upon approach- 
ing the glass transition. The equilibrium tagged-particle 
correlators inherit these properties from the host-liquid 
correlator (j}(t), shown in Fig. [l] as dash-dotted (liquid) 
and dotted (glass) lines. Linear response is the regime for 
small Fcx where the (l>^{t) are still close to their Fox = 
limiting cases, as in this case, C defined through Eq. ^ 
remains force-independent. 

Increasing F^x in the schematic model leads to a de- 
crease of the plateau in the (t>%{t), and to a decrease of 
their a-relaxation time. It can be argued [5] that this cor- 
responds to the fact that the localized probability den- 
sity for the location of the probe particle continuously 
broadens when increasing Fex- For large enough external 
forces, the plateau vanishes completely, corresponding to 
delocalized probe motion. This allows to define a critical 
or threshold force F^^. For glassy states, this indicates 
the force needed to locally melt the glass surrounding the 
probe. In the liquid, the nearest-neighbor cages giving 
rise to glassy dynamics still persist over a time scale t'^, 
and around F^^ these are broken faster and more effec- 
tively by the applied force than by thermal fluctuations. 
For the parameters used in Fig. [T] the following analysis 
confirms F^x ~ 6.5735. 

At still larger forces, the correlator (/)J|(i) in Fig. |lj 
shows oscillatory behavior. This can be connected to 
a finite average probe motion in the delocalized state, 
as we will investigate in more detail below. For (/'^(i), 
no such oscillations are seen, as on average the probe 
will not move perpendicular to the direction of the ap- 
plied force. The oscillations are a clear signature of 
the non-equilibrium nature of the dynamics, since for 
colloidal dynamics, the negative semidefiniteness of the 
time-evolution operator rules them out in equilibrium 

[mini EH- 

Figure [2] presents curves for the real part of the probe- 
particle correlation function corresponding to fluctua- 
tions in the direction of the force, Re(/)j|(t), in a double- 
logarithmic plot. Compared to Fig. [T] a state even closer 
to the glass transition has been taken, in order to bring 
out more clearly the different time scales and the as- 
sociated relaxation laws for the correlators. To explain 
these is the aim of the discussion pursued in this paper. 
An interesting feature brought out by the analysis pre- 
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FIG. 2. Real part of the probe-particle correlation func- 
tions Re0||(t) for the schematic model for force-driven mi- 
crorheology, for a distance to the host-liquid glass transition 
£ = ±10"^ and forces Fex = -F^x(l + 5) with F^^ ^ 6.587 
and 5 — —1, —0.01, and 10~^, from top to bottom. Model 
parameters are v" = 30, k± = 0.5, and fcy = 1 (others as men- 
tioned in the text). Solid (dashed) lines show the results for 
the liquid (glass). Dash-dotted (dotted) lines show the cor- 
responding host-liquid correlator (ji{t) for the liquid (glass). 
Symbols mark the time scales ta and t'^ for the host correla- 
tor (diamonds, cf. Eq. ( |14[ )), tg (squares), ti/2 (circles), and 
tLg (triangles) for the probe correlator (Eqs. (291, (30 1, and 
(32|). 



sented below is the fact that all the probe-particle cor- 
relation functions are, asymptotically close to the glass- 
and the delocalization transition, proportional to each 
other. Hence, a discussion of the time scales seen in 
Re(/)j|(t) suffices. They are marked in Fig. |2| by vari- 
ous symbols, and we will come back to their aiscussion 
below. 

The behavior of the correlation functions shown in 
Figs. [T] and [2] gives rise to a strongly nonlinear signa- 
ture in the friction coefficient. Figure [3| shows the fric- 
tion increment A<^ calculated by Eq. Q as a function 
of external force F^x for several values of e; other pa- 
rameters were chosen as in Fig. [l] Qualitatively, the re- 
sulting C(^cx) agree with those discussed in Refs. [HI HO] 
in conjunction with experimental and simulation data. 
For Fcx —7- 0, a linear-response regime is recovered where 
C depends only weakly on i^ex- C(^ox 0) increases 
strongly with decreasing |£:| from the liquid side. This is 
the manifestation of the glass transition, where the equi- 
librium mobility of the tracer particle vanishes, as long 
as the coupling between probe and host liquid is suffi- 
ciently strong. The possibility of a decoupling of tracer 
motion from the host even in equilibrium (realized, e.g., 
by tracers of sufficiently small size in binary mixtures of 
hard spheres [35]) will not be discussed in this paper. 

As fex is increased, deviations from linear response 
set in quadratically with the force, AC « AC(Fox — >■ 0) — 
aF^^, with a prefactor a depending on e. A steep descent 
of AC is then seen in Fig. [3] around the threshold force. 



FIG. 3. Excess friction experienced by a probe parti- 
cle, as a function of the applied force _Fcx, calculated in the 
schematic MCT model according to Eq. (jlj. Curves from 
bottom to top correspond to distances to the host-liquid glass 



transition e = —1, —0.1, —10 



-10"^ (liquid; sohd lines), 



(solid), 10 0.1, and 1 (glass; dashed lines). Other parame- 
ters are chosen as in Fig. [T] 



Fox ~ ^cx- glass, the point F^^ = F^^ marks the 

divergence of AC as one approaches the critical force from 
above. Fox — >■ F^^ + 0. At forces much larger than Fox, 
a second plateau is observed in AC. Intuitive reasoning 
might suggest C(-fex — > oo) — )■ 1 (the solvent friction in 
our units), hence AC = in this window. This is not seen 
in experiment or simulation. It was argued in Ref. |10j 
that within ITT-MCT, the fact that AC > for Fox oo 
can be understood by accounting for contributions to AC 
stemming from a =_L, i.e., fluctuations in the direction 
perpendicular to the force. Note that no reference to 
a suspending liquid is made in this argument, although 
hydrodynamic interactions mediated through the solvent 
are expected to greatly influence the friction coefficient 
for large forces in real colloidal suspensions. The ratio of 
the large-force plateau to the linear-response value has 
been determined in a low-density expansion by Brady 
and coworkers [531 [M] to be AC(Fox oo)/AC(Fox ^ 
0) = 2; this is confirmed by simulations of dilute host 
liquids. In the schematic model, the parameter fi serves 
to reproduce this ratio when one assumes the coupling 
coefficients (wi,W2) and (wi,!"!) to approach zero at low 
densities: dropping all memory kernels from Eqs. (pi) and 
([3]), the correlation functions become 4'{t) = (f>iXt) = 
exp{—t) andRe0||(f) = exp(— i) cos(Fexi). InEq. ^ this 
yields ACo = (4 + fiF^J/{8 + 2F^J, which interpolates 
between ACo(Fox 0) = 1/2 and ACo(Fox ^ 00) = fj,/2. 
In Fig. [3] this is not seen, since we kept v'^ = 30 fixed. 
While fits to data and the microscopic theory suggest to 
change v" as a function of density [l^, for a discussion 
of the features close to the glass transition density we 
assume w D*'(e = 0) for simplicity. 

We now focus on the analytic discussion of the above 
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observed crossovers. For this purpose, let us employ a 
more compact notation of Eq. (|3| . We introduce a vector 
of probe correlation functions, (p^it), with components 

(0f(i),0^(i),0^W)^ := (Re^|(i);im^|(i),01(O)^, and 
initial conditions 0^ = (1,0,1)^. We will refer to the 
components of this vector with Latin indices j S {1,2,3} 
(recall that Greek indices, a E {||,_L}, label the direc- 
tions). Let us further introduce a matrix characterizing 
the short-time motion, 



U! :- 




(5a) 



and a matrix constructing the memory kernel vector. 



Vn 
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(5b) 

with = 1 + (naFcx)'^- The complex-number mul- 

tiplication appearing in the convolution integral is then 
represented in our matrix notation by a bilinear (but non- 
symmetric) mapping C|[x, y] :— ■ A4^^ C . ■ y, with 



a. 



a. 



It will occasionally be useful to introduce the sym- 
metrized version of this mapping, I?'* [x, y] = Cf [x, y] + 
C^[y,x\. Applying the Laplace transform, 4'{z) = 
i ^ e'^^*(f){t)dt, the equations of motion Eqs. (2a) and 



(3a 



are rewritten as 

z4){z) - [l + z(j){z)){iz + zm{z)) ^ Q (6a) 
for the bath correlator, and 



-1 (z0'(z)+0^) +Z0'(2) 



zC 



for the probe correlator. 



\z),z'^ (z) + 0^ -0 (6b) 



III. LONG-TIME BEHAVIOR 

We assume the long-time limits of all correlation func- 
tions of the schematic model to exist. In view of standard 
features of structural relaxation dynamics close to a glass 
transition, this is not unreasonable. Then, the Abelian 
theorem |25j. 



leads from Eqs. ^ to 

fl{l-f)=V,f + V2f 

for the bath correlator and 



(8) 



Q^A'-f~fc![f_\n, (9) 

for the probe correlators. Here we have defined 

g-r--^fcv_^^]-r. (loa) 

Explicit evaluation leads to 
A'= U||K||Fexi;f/ -i\\vU-l i\\n\\F,^v-J \ . (10b) 







Equation ([s]) describes the well-known bifurcation sce- 
nario of the glass transition within MCT [15]: it is an 
implicit nonlinear equation for the nonergodicity param- 
eter /. Possibly many different solutions of this equation 
exist for general models, with zero being an obvious one. 
It can be proven [13] that within equilibrium MCT, / 
is always determined by the largest positive real solu- 
tion. Since for ui,2 — > 0, only f — survives, and for 
vi^2 — >■ oo, / > holds, the MCT solutions changes at 
some bifurcation point. These points define a hypersur- 
face (wi,c,i^2,c) through vi,c = 2y^t^ — V2,c (restricting 
to 1 < 1)2, c £ 4). The bifurcation is identified as the ideal 
glass transition. Generically, in Eq. ([8|, a jump in / oc- 
curs; this is called a type B transition in the literature. 
In other words, /c > holds for the critical nonergodicity 
parameter evaluated right at the transition. Mathemat- 
ically, one deals for generic mode-coupling models with 
bifurcations of the class Ae [13] according to the Arnol'd 
classification [26^; in the present case we are only con- 
cerned with the nondegenerate A2 bifurcations displayed 
by the F12 model. The bifurcation points are identified 
by recognizing that, as two branches of solutions of the 
implicit equation for / coalesce, the equation is no longer 
invertible (and the implicit-function theorem violated). 

Equation ([9| describes localization respectively delo- 
calization of probe particles: if / = 0, also = 0, i.e. 
in a liquid, all probe particles are delocalized (and able 
to undergo long-range motion, as required for a liquid). 
If / > 0, sufficiently large {vf,V2) will for F^^ = lead 
to > 0, so that a probe is localized in the glass if it 
couples strongly enough. As noted above, we will not 
discuss the weak-coupling limit where the probe remains 
delocalized in the glass even for Fox = 0. 

For non-zero Fex, a further bifurcation point will de- 
scribe probe delocalization in the glass through external 
force. It can be found by demanding that the implicit- 
function theorem is violated for Eq. ([9]), det A" = 0, lead- 
ing to a biquadratic equation. 



/ := hm cj){t) 

t^OO 

r :- hm cf>'{t) 



■ lim z(j){z) , 

2— >+iO 



(7a) 



- lim z(j) (z) , (7b) 

z^+iO — 



^ 44(^eJ'-(«?/ - 1) (4 K/ + 1) + 4) (^ex)' 
+ {vtf+l)({vtf-lf-iv^2ff) ■ (11) 
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FIG. 4. State diagram for the probe particle in the schematic 
model. Upper panel: critical force F^x as a function of dis- 
tance to the host-liquid glass transition e, for fixed as in 
Fig. [l] Lower panel: critical force as a function of probe- 
to-host coupling strength v" f . Solid lines mark the bound- 
ary between localized and delocalized probe states; a vertical 
dotted line in the upper panel indicates the host-liquid glass 
transition. Dashed lines are f^x obtained for the model with- 
out taking into account the anisotropy of the probe-particle 
fluctuations, V2 = 0. 



Solving Eq. (11) to obtain F^^^ as a function of u"/, 



we obtain a decomposition of the parameter space for 
our schematic model. The result for exemplary parame- 
ters is shown in Fig. [4j The lower panel shows F^^ as a 
function of w*/ directly. Typically, f will change im- 
plicitly as the glass form factor / changes with changing 
e; this variation is shown in the upper panel of the fig- 
ure for the parameters chosen as above. In Fig. |4j also 
the result after setting u| = (as implicit in the model 
of Ref. 0) is shown for comparison (dashed lines). For 
this aspect of the discussion, both models are very sim- 
ilar. For v'^ f < 1, no real solution of Eq. (11) exists, 



and the probe particle is always delocalized. This is triv- 
ially the case in the liquid, where / = 0. In the glass, 
/ > holds, and by assumption we restrict ourselves 




FIG. 5. Solutions of the implicit set of equations Eq. (jof, 
/i = Re /I (sohd lines), /| = Im/|| (dash-dotted), and 
/I = f1 (dashed), as a function of the reduced applied force, 
Fcx/i^"/)- Parameters are = 30, / = 1 — l/\/2 corre- 
sponding to the glass-transition point of the F12 model at 
V2 — 2. The solutions identifled as the long-time limits of the 
schematic model, Eq. (|3|, are shown as thick lines. The inset 
shows the full set of solution branches, including one for /| 
that is cut off in the full figure for clarity. 



to the case where the probe particle also becomes local- 
ized without any external force, so that the jump from 
zero to fc is big enough to render fc > 1 at the glass 
transition. As exemplified in the upper panel of Fig. [4j 
this leads to a nonvanishing critical force F^^ > for the 
delocalization of the probe at the glass transition, iden- 
tifying the regime Fox < Fgx ^ the one where the probe 
remains localized inside the glass even under the action 
of an external force. Beyond that, / increases with e, 
asymptotically as / — /c ~ \/e. Thus also F^^ increases 
in the glass. A full analysis of Eq. (11) shows that a 



further line of solutions is present inside this localized 
regime. Based on the following discussion of f^, we do 
not assign physical significance to this. 

Solving Eq. ^ numerically, we obtain a set of solu- 
tions to this implicit equation. These are shown in Fig. [5] 
where we have highlighted those in bold that we will use 
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in the following as the physical ones. For Fox < they 
are the continuation of the known solution /^(Fex = 0). 
In particular, /| — > as — > since the imaginary part 
of the correlation function has to vanish in equilibrium, 
and fl = /I for Fcx — 0. The physical branch crosses 
zero at F^^, and for larger Fcx, no solution emerges for 
which both real parts, ff and /|, are non- negative real, 
except = 0. We hence are dealing with a continu- 
ous probe-delocalization transition: the probe-correlator 
nonergodicity parameters /" do not exhibit a jump at 
Fg^. Note that the choice of solution branches is non- 
trivial: in standard equilibrium MCT, it can be proven 
that the solution that is the largest nonnegative for all 
components has to be chosen [^. The proof relies on 
mathematical properties of Eq. ([8| and its proper gener- 
alizations, that are not ensured in general in Eq. The 
number of solution branches to Eq. ^ may depend on 
the parameters, and in particular the choice of the vf and 
^2 ■ Fig. |5] shows a case that is typical for the description 
of microrheology data. 

It can be seen that at F^x = F^^, the largest real eigen- 
value of the matrix crosses zero, and that this eigen- 
value is always non-degenerate for the restriction to the 
parameters introduced above. Thus the critical force is 
determined by a bifurcation in Eq. Q with a single criti- 
cal direction (co-dimension = 1^ It is in this sense 
similar to the known continuous delocalization transi- 
tions discussed within tagged-particle models of MCT 
|15j . These are quite different from the standard MCT 
bifurcation, as we are not dealing with the coalescence 
of two solution branches, but with the crossing of one 
particular solution with the /J = branch. F^^ is well- 
defined since setting any of the fj = in Eq. ([9| also 
demands that the other f^_^j = 0. 

For the later discussion we define reduced distances to 
the critical points of the model, e and S. We set 



V = {vi,V2V = V, + e-b^, (12a) 
W=ivlv'2,F^x)'^ = W'' + 6-b''\ (12b) 



The vectors and If''^ can in prinicple be chosen arbi- 
trarily. In particular, the weak-coupling limit of small v'^ 
is contained in the following derivation. For the sake of 
simplicity, we will choose for all explicit calculations, in 
agreement with above. 



he = ivi_c,V2,cV 



(12c) 



With this choice, 5 becomes a reduced force, such that 
(5 < corresponds to Fcx < F^^, i.e., the localized regime, 
and (5 > to the delocalized regime. Recall that sub- 
scripts c refer to quantities evaluated at e = 0; we reserve 
superscripts c for quantities evaluated at S = 0. 



IV. THE BETA-SCALING LAW 

A. Scaling Functions 

Having identified the long-time limits of the correla- 
tion functions in the glassy-host, localized-probe regime, 
we now turn to a discussion of the dynamical correla- 
tion functions ^(<) and (t)%{t). We anticipate that upon 
approaching a critical point, e — for small 5, the cor- 
relators will stay close to an intermediate plateau given 
by / over a time window that increases as the distance 
parameters approach zero. Thus, (t){t) — / and (t>'^{t) — f 
can be identified as small parameters on a dynamical time 
scale that will be determined below. Here, / is a param- 
eter introduced in Ref. [27 , approximating / for e > 
and continuing the latter smoothly to e < such that 
lim£_^±o / = fc- For the F12 model, where we only con- 
sider A2 bifurcations, / is given by the unique positive 
solution of vi + 2v2f ~ 1/(1 — /)^ — that obeys f = fc 
at the critical point. We define / and a corresponding 

matrix A by Eqs. ^ and ( 10 1 with / replaced by /. 

For the host correlator (jj given by Eqs. ([2]), the asymp- 
totic analysis has been worked out in detail. Splitting 
(f){t) = f + G{t), there holds an expansion in terms of 
the small parameter a: G{t) — c„hg„{t) + 0{c^) with 
c<T — \/M = ©(-y/jej). The leading order correction to 
the plateau is called /3-correlator gcr{t). There appears a 
critical amplitude h, connected to the critical eigenvector 
arising in the bifurcation analysis of Eq. ([s]) - in the one- 
correlator schematic model, this is just a prefactor, set to 
h — (1 — fc) by convention. The function g{t) describes 
relaxation to and from the plateau (visible e.g. in Fig. [T]) 
as asymptotic power laws. 



g^{t) cx (tAo)^" 
g^{t) cx -{t/uf 



for <o < i < ^<T 
for -^t-^C 



(13a) 
(13b) 



where the latter, called the von Schweidler law, only oc- 
curs inside the liquid, e < 0. Here, diverging time scales 
have been identified, 



tn 



tJlcr 



l/(2a) 



t' 



tn a 



l/(2a) + l/(2fc) 



(14) 



where cr is a distance parameter obeying cr ~ e as e — > 
0. For the schematic model and our definition of e, one 
gets a = efil + /=)(1 - f'') « 0.38e. The microscopic 
time scale to is fixed by the short-time motion of the 
correlation function. The power-law exponents a,b > 
are solutions of 



r(i-a)2 _ r(i + 6)2 

r(l - 2a) ^ r(l + 26) 



A. 



(15) 



with the exponent parameter A given by A = 1/ y/v2 on 
the bifurcation manifold ^15, . 

Introducing a rescaled time, t = t/t^, the /3 correlator 
is found to obey the /3-scaling equation written with the 
scaled Laplace frequency z — zt^ as 



Tl + A%|(z) + (zg±(z))2 



0, 



(16) 
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where g±{t) are the scahng solutions for cr ^ that de- 
pend only on the sign of a. They obey g± (t) ^ for 
t — > 0, and g-{i) ^ —Bt^ for i — > oo. The asymptotic 
form of the relaxation to and from the plateau thus is 
given by a scaling law with only A as a non-universal pa- 
rameter. The positive constant B ~ 0(1) depends on A 
and has been tabulated For e > 0, g+{t) approaches 
a constant as t — > oo, giving the correction to the critical 
plateau value. 

In the same spirit we now write 



(17) 



To arrive at an equation for g^ that is asymptotically 
valid for cr ~ £ — > 0, we follow the standard procedure 
laid out in Ref. |15) : introducing z = zt^^ we collect in 
Eq. ( 6b ) the leading terms in Ca , making use of Eq. ^. 
We arrive at 

= {A'-zgl{z)-fmi\S9lim 

+zg.m'[r,^i-r]} 

+ Ca [Qf[zgagl{z),(l>_l-T] + Zg^{z)Cf[~f Jgjz)] 

+ /Cl5g:(5),i5:(5)]} +0(c^,M.)-^) (18) 

Note that since the exponent a < 1/2 [15], the dropped 
terms (c„t„)^^ ^ |(j|i/(2a)-i/2 g^j-g indeed of higher order. 

The leading order for e 
linear equation system, 



in Eq. (18) results in a 



^:-£l(i)=5± it) 11, 



(19) 



where the linear mapping is defined through ~ 
4^ ~ /c2^c[/c''] ^^"^ inhomogeneity reads — 
— C'*[/^,(/)^ — /^] . For state points far enough from the 



critical force, Eq. ( 19 1 determines 5^ {i) by the host liq- 
uid g±{t). This expresses that in such cases, the asymp- 
totic dynamics of the tracer is governed by that of the 
host liquid, i.e., 



9lii) = hlg±{i), 5^0, 



with an amplitude given by 



(20a) 



(20b) 



In Fig. [6] we show exemplary results for the corre- 
lators <t>^{t) close to the glass transition, in a double- 
logarithmic plot of \(l)^{t) ~ to exhibit the asymptotic 
power laws. Two states corresponding to small \e\ were 
chosen that exhibit a large window of validity for the 
asymptotic Eq. ([l8|. As dash-dotted and dotted hnes, 
the corresponding (/)(i) are shown for e < and £ > 0, 
respectively. They follow the asymptotic critical power 
law, (j){t) — ^ h{t/to)~°' over several decades in time; a 
numerical estimate yields to ss 0.146 for the parameters 
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FIG. 6. Correlation functions </>''(f) close to the glass transi- 
tion point, for various external forces, plotted as the leading- 
order deviation from the plateau value, |G*(t)| = |0^(i) — .f"\ 
for e — > 0, as a function of t/ta- Parameters as in Fig.|3] solid 
lines correspond to liquid states, e = —2.5 x lO"'^", dashed 
lines to glassy states, e = 2.5 x 10"^". Various forces Fox 
are shown with 5 = (Fcx - Fox)/f?x = 10"", n = 0, 1, 2, 3, 
4, and 5 as labeled. (For the glass, only n = and n — 1 
are shown.) Dash-dotted and dotted lines indicate the corre- 
sponding host /^-correlator jG(t)|, for the liquid respectively 
the glass, in the panels showing the real parts. For the imag- 
inary part, dotted lines indicate the critical law of the host 
correlator, G{t) ~ /i(t/to)"", and Kj/^ ~ 3.985/A times this 
critical law. 



chosen here. The £ < curve also exhibits von Schwei- 
dler's law, (f){t) — ^ ~Bt^ at large rescaled times 
i — t/t„. With the distance parameter chosen here, we 
get (J — 6.7x 10"^^ and consequently = 8.07x 10^^; the 
parameter B is tabulated [3S] and estimated as -B = 0.68. 

For forces far from the critical threshold, e.g., 5 = —1, 
Eq. (20) is valid. This is demonstrated by the n — 



curves in Fig. [6] that correspond to the force-free case: 
the probe-particle correlation functions closely follow- 
ing the host-liquid curves over all time, up to a fixed 
amplitude. Evaluation of Eq. (20 1 in this case yields 



hi (0.39, 0, 0.39)^, which is easily verified in the figure. 
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In this sense, the tagged-particle dynamics is governed 
by the dynamics of the host hquid. This well known fact 
is the basis for various approximations relating single- 
particle motion and collective dynamics close to the glass 
transition. For n = 1, we get « (3.99,2.83,2.54)^, 
and this is also verified by Fig. [6) 

Approaching F^^^ shown by curves with S = — 10~" 
closer to zero (larger n) in Fig. [6j the coupling found for 
6 = breaks down. Most obviously, in the regime of the 
von Schweidler law, t ^ 1, the curves for n ^ decay to 
zero more rapidly than the host correlator, indicated by 
a plateau in Fig. |6] (dashed lines). This regime will be 
discussed later. 

More subtly, for t ^ 1, the probe correlators for 
^ex 7^ still appear proportional to the host liquid 
correlator, but with a different prefactor than the one 
given in Eq. pa. For n = 5, Eq. |20| yields w 
(4.16, 3.98, 2. 05y7 not compatible with the asymptotes 
seen in Fig. [6j We address this now, discussing the limit 
(5 -> of Eq. ([18]). 

At the critical force, S = 0, we have detA'^''^ = and 
concomitantly lim5_j.o = /^''^ = 0, as discussed above. 
In that case, also det I/^''^ = 0, rendering Eq. ( 19 ) singular 
and Eq. (20) invalid. In leading order in S, 5* (t) must be 



in the kernel of det A^''^. In our model, this matrix is only 
simply degenerate, and we denote by If''' a corresponding 
null-eigenvector. Keeping the ansatz that g^(t) depends 
linearly on gait), we modify Eq. ( |20a ). 



9:it)=ll''gaii) + 0{ca,S), 



0. 



(21) 



Differentiating Eq. (|9| with respect to 6 at S = 0, again 
using that / ' = 0, we find A ' {d/dS)f = and hence 

{d/dS)f ' = -qh''" with some constant 77. Here and in 
the following, expressions like d/d6 f^''^ are to be read as 
lim5_>o d/dS f". 

Accounting for the next-to-leading order terms in 
Eq. (18) and using /J''^ = 0, we get from taking the 
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n = 5, 




derivative with respect to 6 



FIG. 7. Normalized ratio of probe-correlation functions to 
that of the host liquid, X,(t) = [(?!.,'(t)-//]/[<^(t)-/]/(u;^Ma), 
cf. Eq. ( |25| j as a function of f = t/tc- Parameters are chosen 
as in Fig. IgI only liquid curves, e = —2.5 x 10"^" are shown, 
for 5 = —10"" with n = 0, 1, 2, 3, 4, and 5 (indicated by the 
labels). The dotted line indicates A^/As = A = \/2- Squares 
in the upper panel indicate the time scale ts- 



= {zga{z))S X 



dA 



dS 



+ ca {c^■^[A^^^0^]z^(z)-h/C^■^[/^^^^A^^=](z5.(^))2} 

+ 0{cl,{cata)-\6^cJ) (22) 

In the following, we will suppress the indication of omit- 
ted higher-order terms as far as they are already denoted 
here. The term containing {d/dd)A can be eliminated 
in favor of Cf''^ by looking at the Taylor expansion of 
Eq. (|9|, which yields 



= A 



= 2 d(52 



dA 
^~d5~ 



(23) 



Multiplying with the left null-eigenvector h ot A , the 
first term vanishes, and inserting into Eq. (22) provides 



+ft'C''^[h^'',h^'^]{zgAz)f} (24) 



Let us discuss the relevant limits of Eq. (24). We as- 



sume that C^'''[h''''',If'''] 0, which is generically the 
case. Note that in writing Eq. ( [2T| ) , the normalization of 
the eigenvector hf''^ is left undetermined. In the limit 
£ — >■ 0, this normalization is now fixed by requiring 
the first bracket in Eq. (24) to vanish. If we introduce 
through hf''^ = Au the normalized eigenvector u and the 
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amplitude A 



, we obtain 



(25) 



where the subscript recalls that this result is valid for 
taking the limit (5 — > after e — > 0. The amplitude 
As does not depend on the actual path along which the 
hypersurface of delocalization transitions is crossed. 

Taking (5 —> first, wc require the second bracket in 
Eq. (24 1 to vanish. Letting e — after that, we can 



further make use of Eq. ( 16 1 to arrive at 



A, = 



Xfch 



A^ 
A 



(26) 



Equations (25 1 and (261 together with Eq. (21) imply 



that one recovers both the critical power law, 5^(i) ^ 
t"", and the von Schweidler law, 5* (t) ~ —i^, in the dou- 
ble limit £ — > and 5 — 0. The ratio of the correlation 
functions, however, depends on the order in which this 
double limit is taken. To exemplify this, define the ratio 
X^{t) = {r,{t)-f^)/mt)-f]Asu,). Results for typical 
parameters are shown in Fig. [7j According to Eqs. (25) 
and ( 26 ) , the ratio approaches unity for i; — >■ 00 at not too 
small 5. This is demonstrated by the n = curve in the 
figure. Lowering S ~> 0, one notices a second plateau in 
Xj{t) = 1/A that becomes more pronounced for smaller 
\S\. For Fig. [rj the result is only valid for t < ta- The 
ranges of validity of the various scaling predictions will 
be discussed below. 



Another scaling result emerges if instead of Eq. (21) 



we consider the case that ga{t) is of higher order and can 
be set to zero in Eq. (18). We will identify a time scale 
where this is admissible below. Requiring the leading 



order in Eq. ( 18 1 to vanish again leads to the requirement 
cx hf^''^^ the null-eigenvector of A ' . The expansion 



Eq. ( 22 ) reduces to 



d6 




(27) 

The limit e — > is not meaningful in this equation, since 
the first bracket then only admits the trivial solution. 
Letting, however, 5 — )■ 0, we gain C'*''^[zg^(z), zg^(z)] — 
0, which is consistent if we set 



' (f) = cr^/^hf^'^ . 



(28) 



The constant C has to be determined by matching the 
various asymptotic expansions. 



B. Time Scales 

The different results derived above present the asymp- 
totic behavior of the probe-particle correlation functions 
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FIG. 8. Correlation functions <j)"{t) — f" multiplied by -v/t, 
for e = 10~®, and other parameters as in Fig. [g] Reduced 
forces are S — — 10~" with n = 0, . . . 6. The dash-dotted line 
shows the corresponding host-liquid correlator. 



in various limits. These correspond to different time 
scales for which the results hold, which we discuss now. 
Recall the Laplace transform of a power law g{t) = 
yields zg{z) = —V{\ — x){—izY. For short rescaled times, 
t — > 0, the host-liquid /3 correlator assumes the form 
gcr{i) ^ t^°. In Eq. ( |24| ), we then identify a time scale 
ts that separates the two scaling limits discussed in con- 
nection with Eqs. (25) and (26 1: for t ^ i^, the second 



curly bracket in Eq. (24) dominates, and we recover the 
(5 — > regime. For t ^ tg, the first curly bracket domi- 
nates, corresponding to the £ — >■ regime. Balancing the 
power-law exponents, we get is — 0{l)\a/S'^\^/'^^'^\ or 



(29) 



Considering the case where ts S> to-, we have to dis- 
tinguish the two signs of £. In the glass, £ > 0, ^o- at- 
tains a long time limit of 0{\/a) so that we may consider 
Eq. (27l to be vahd on a time scale ii/2 = 0{l)\a/S^\, or 



(30) 
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For times t ^ ti/2, the trivial solution of Eq. ( [27| is 
obtained. Hence the probe /3 correlators G^{t) show a 
nonalgebraic decay to zero for times t ^ ii/2- Matching 



This yields the 



this solution, Eq. (28), to the asymptote obtained above, 
we find that the constant C obeys C = C" • |cr|i/2-i/(4a)^ 
where C" is independent of a and 6. 
scaling law 

G^(^)(x|cTrl/(4'^)|<5|(^Al/2)-l/^ 



< t < ■ (31) 

We verify this scaling law in Fig.|8[ for e > much larger 
than considered in Fig. [Hj G^(t)\/t is seen to approach a 
constant at times long compared to the ultimately expo- 
nential relaxation jT3j of G{t). 

In the liquid, £ < 0, we can balance terms of 0{S) 
with those of 0{c„) inserting the von Schweidler law for 
the host-liquid (3 correlator, ga{t) ^ — C'(l)i~^, provided 
that i 3> 1. As a result, this asymptotic power law can 
be seen in the probe correlator on a time scale 



(32) 



Taken together, Eqs. (14), (29), (30), and (32) define 



five time scales determining the asymptotic behavior of 
the f3 correlation functions for the probe pulled by an 
external force. Two of them are inherited from the glass- 
transition dynamics of the host liquid; tu determines the 
time scale for the relaxation around the plateau, while 
t'^ sets the time scale for the final decay in the liquid. 
These time scales are marked as diamonds in Fig. [2j we 
have seen that it suffices to study the correlator Re (/<| (t) 
in the asymptotic regime, since the other (t>j(t) are con- 
nected via the critical amplitudes /i*''^ to the same time- 
dependent laws. 

The time scale ts (squares in Fig. [2]) separates the 
regime that is dominated by the proximity to the glass 
transition (cr — > 0) from the one dominated by the prox- 
imity to the delocalization transition ((5 — )■ 0) . The curve 
for 6 = —0.01 exemplifies the case where ts <C so 
that the dynamics for t ^ ts is dominated by the host- 
liquid glass transition. In this case, the probe-particle /3 
correlator is proportional to the host-liquid one, with an 
amplitude that changes by a factor 1/A at t « ts. For 
t ^ ts, Eq. (26) holds for this amplitude, while a,t t ^ ts 
it drops to the value determined by Eq. ( 25 ) . This ex- 



plains the behavior of the correlator ratio Xit) shown in 
Fig. [7] There, squares mark ts; for the particular cases 
shown, ts < t„ holds, so that the ratio for t > t„ \s 



given by Eq. (25), viz. Xj{t) — 1, as long as the probe 
correlators do not decay to zero. 

At times t ^ t^r ts, the probe-particle corre- 
lator exhibits force-induced decay that couples to the 
von Schweidler law of the host liquid: this is exemplified 
for 5 = —0.01 in Fig. [2] for i w ^ (marked by a trian- 
gle). There holds t'^ ^ (x. \S\^^^ , so that approaching the 
delocalization threshold, (5^-0, this force-induced decay 
is accelerated according to a non-trivial power law. 



The curve for 5 = 10~^ in Fig. [2] on the other hand 
exhibits the case ts ^ ta, and the dynamics for t <^ ts 
is dominated by the probe-delocalization limit. In this 
case, the time scale becomes relevant (marked by 
a circle in the figure), and the probe correlator decays 
according to a power law. 

Both ts and ti/2 diverge as the force threshold is ap- 
proached, (5 — > 0, while the host-liquid time scales t^ and 
t'^, as well as the time scale for which the probe motion 
couples to the host's von Schweidler law, t'^ g, diverge for 
£ — ^ 0, i.e., upon approaching the glass transition. From 
the power laws, one infers the relevant scaling combina- 
tion: for |cr/5^| 1, the dynamics is determined by the 
force-induced delocalization, while for |cr/(5^| ^ 1, the 
glass-transition dynamics describes the relaxation. 



V. THE a-SCALING LAW 

In the liquid, the decay of the host correlation function 
from the plateau towards zero can be discussed on the 
time scale t'^ by another scaling law, known as the en- 
sealing law. There, a master function can be derived for 
the limit cr —>■ in the form [T5] 



0(i) X F(VC) . 



(33) 



The function F is independent of a. As a result, the 
final relaxation of the correlators displays a superposition 
principle, allowing the long-time part of the correlators 
to be scaled to a master curve by scaling time with t'^. 
We derive a similar law for the probe correlators ^*(i). 

More formally, we introduce cr-independent scaling cor- 
relators F and F^ by 



F(t) :== lim (j)(t) , 



Flit) := lim 



^%t), (34) 



with reduced times t = t/t'^ and t = t/t'^g. From 
Eqs. we obtain after taking the scaling limit £ — )■ 
with t^,t'^ g 00: 



F{-z)=T^[F]i-z)-il + -zF{-z)), 
11(2)= CI \FF!iS),SFl{2)+^l 



(35a) 
(35b) 



Here, J^[F] denotes the memory kernel of the host-liquid 
model, evaluated with the scaling correlator. The initial 
conditions for Eqs. (35) are obtained by matching them 



to the long-time limit of the /3-correlation regime in the 
liquid, i.e. the von Schweidler law. Hence, F(0) = /c, and 
limf_j.o(F(t) — fc)t~^ = —B. This leads to the standard 
a-scaling regime discussed in detail earlier [15j . Likewise, 
we identify F^(0) = and lim.^^o (F^ (i) - Qi-' = 
—B'^ where is trivially related to B via the critical 
amplitudes. As the plateau value vanishes linearly 
with 5 upon approaching the delocafization transition, 
it is reasonable to introduce the probe-particle a-scaling 
correlator 



'(t):= hm m)m. 

o— )-0 



(36a) 
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to which corresponds the a-scahng law 

Um nm0^(t)/|<5Hr(i/4,) 

(5->-0- e->0 — 



(36b) 



It remains to be shown that F_^{i) indeed is indepen- 
dent on cr and S. To see this, expand Eq. (35b), 



= 4"'" = 



SF {z) 



dS 



sfF%S),ll\}+gf-0{S). 



(37) 



This suggests the ansatz F^^i) = F'{i)h'f + 0{5), and 



multiplying Eq. (37) with the left-null-eigenvector 



yields an equation that is independent on 5: 



= -SF'{S)h^ 



dA' 



d6 



U~zF^{~z))%T -CTlhTAT] 

-B^~zfF^{~z)hT -Crm^ 



(38) 



Consequently, F_^ is invariant under rescaling of both 
variables 5 and e. It depends, however, on the direction 
of the path crossing the delocalization transition with 
respect to the transition hypersurface. 

The scaling property of the probe correlators as a func- 
tion of (5 0^ in the limit e is exhibited by Fig. |9] 
Here, some very small e < has been chosen, and the 
correlation functions of the probe particle are plotted as 
4>a{i)/\^\ foi' various 5 < Q. Reducing \5\, the a-scaling 
regime is entered, as seen for the n — 2 through n = 5 
curves: as functions of rescaled time t = t/t'aS' 
correlators agree with the scaling function F^{t) at long 
times. The latter shows an initial plateau, followed by 
a relaxation on a time scale i — 0(1). Qualitatively, 
this scaling is similar to the usual a-relaxation scaling 
law one finds without external force upon varying e — >■ 0. 
There, however, the discontinuous nature of the MCT 
transition predicts that asymptotically, the plateau seen 
in the a correlator is constant, while for the delocaliza- 
tion transition approached as J — )■ 0, the plateau scales 
with the distance |^| to the transition point. For n = 
and n — 1, the distance \6\ to the delocalization transi- 
tion is too large, and the scaling prediction of Eq. ( 36 ) no 
longer holds, 
curve for n = 

^2 , 



as it requires both 5 — > and e — >■ 0. The 
= 6 also shown in Fig. [9] corresponds to the 
regime <C 1, and no longer satisfies the require- 

ment e in Eq. ( 36 ) ; it hence deviates again from the 
scaling law. 



VI. THE FRICTION COEFFICIENT 

A. Above the Delocalization Threshold 

We now focus on a discussion of the friction coefficient 
AC in the asymptotic regime. Let us start by considering 




-12 -8 -4 

logio t/t'cb 

FIG. 9. Test of the a-sclaing law for the probe-particle 
correlation functions: the (j)" {t)/\5\ are shown as functions of 
rescaled time t = t/t'^ g. Parameters are chosen as in Fig. [T] 
with £ = —10"^^. Reduced forces are given hy 5 = — 10~", 
and n = 0, 1, 2, 3, 4, 5, and 6 as labeled. 



5 > 0, i.e., a force exceeding the delocalization threshold, 
so that even in the glass, AC < oo. 

Figure [TO] shows the friction coefficients from Fig. [3] in 
a double-logarithmic representation as a function of the 
distance to the delocalization threshold, 5 > 0. At large 
(5, the qualitative features discussed before in conjunction 
with Fig. [3] are recovered: a plateau as (5 — > cxd, and an 
intermediate decay for S 1 that we will identify with a 
power law. For S <^ this power law crosses over: for 

liquid states, the curves approach a constant, since the 
friction at the critical force {S = 0) is finite there. For 
glassy states, the friction diverges as 5 — )■ 0+, but with a 
different power law than the one observed for 5 3> 

Consider first the regime dominated by the S ^ 
limit: here, Eq. (28) describes the power-law decay of 
Gf{t) - for cr > 0, valid in a time window 



AC 



< t < 
^ r*i/2 



^1/2- 

C"\c 



Inserting into Eq. (|4|, one obtains 

|l/2-l/(4a) . (^^^^^1/2 . ^ ^^.^^ ^^^^^^^ 
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FIG. 10. Friction coefficient increment Ai^ = C(i^cx) — Co for 
a probe particle pulled with external force Fex, as a func- 
tion of the distance to the delocalization threshold, S = 
(Fcx — F^x)/Fcx, for S > and various distances to the glass 
transition e = ±10"*^ with fc = 5, 6, 7, and 8 as labeled. 
Dashed lines are in the glass (e > 0), solid lines in the liquid 
(e < 0). The curve at the glass transition, e = 0, is shown as 
a dash-dotted line. Inset: curves for e 7^ 0, as a scaling plot, 
AC = AC/lel^'"'-'''*"'"' versus 5 = 5/\e\'^\ 



AC((5) cx |(7|'-'/('°) -,5-1, a > 0, |(5V(t| < 1. (39a) 

In the liquid state, the t^^/^ power law does not appear. 
For 5 ^ we get the (5-independent result in the 

liquid, 

Add) CX |^|V2-l/(2a) ^ ^ ^ 1^2/^1 ^ ^ (39^^) 

Equations ( |39p e xplain the qualitative behavior seen for 
(5 ^ in Fig:]lO| 

For ts ^ tcr, the probe /3 correlator G''(t) for times 
to t ^ tcr follows the power law (io/t)°. Integration 
of Eq. Q leads to 



AC CX r , \6ya\^l, 



(40) 



where the proportionality constant is independent of a. 



This power law is seen in Fig. 10 in the intermediate-<5 
regime, where 6 ^ but still small. 

Having identified these power laws, a scaling predic- 
tion is obtained for the friction coefficient, for forces 
just above the threshold. Introducing S — S/^/\e\ and 
AC = A(^/\£\^^^~^^^'^°'\ we obtain two master curves, one 
for e > 0, and one for e < 0. This is shown in the 
inset of Fig. [TO] The scaling curve for the glass shows 
two power laws, AC ^ S"^ for (5 ^ 1, and AC ^ (5^^^^° 
for 6^1. For our model, a ~ 1/3, so that the ex- 
ponent governing the large- j decay in AC is close to 2, 
as seen by inspection of Fig. [TO] One has to note that 
the identification of these power laws in experiment or 




FIG. 11. Friction coefficient increment AC as a function of 
the distance to the delocalization threshold S, for S < (i.e., 
-fox < J"cx) in the liquid. The distance to the glass transition 
is chosen as e = — 10"*°, with fc = 4, 5, 6, 7, and 8 as labeled. 
Inset: scaling plot AC versus 5 for the same data and also 
including k = 9 and = 10. 



simulation is hampered by the fact that one has to ap- 
proach both £ — ^ and 6 ~> 0. For finite distances to 
the transition, preasymptotic corrections quickly become 
dominant. Fits using the schematic model we discuss 
here have to be performed with varying coupling coeffi- 
cient [TCT. Compared to the case displayed by Fig. 10 
this introduces a regular shift of the curves [29]. 



B. Below the Delocalization Threshold 

We now discuss the case (5 < 0. Here, the probe re- 
mains localized within glassy states, so that for (5 < 
and fj > 0, AC = 00 holds trivially. For ct < 0, the fi- 
nal a relaxation will give the dominant contribution to 
the friction coefficient, which is obtained by inserting the 
a-scaling law into Eq. Q, 



AC«(1-m)/ Fit/QFl,{t/t',)dt 



fi I F{t/t'„)FlS/Cs)dt. 
'0 



(41) 



This expression yields a factorization of AC of the form 
AC - |CT|-i/(2a)-i/(2fc)^(-^) fQj. ^ 0, where Z is uni- 
versal. If 5 ^ i<j, i-c, 1(5^/(71 ^ 1, and additionally 
\S\ ^ 1, we further use that F^s{t/t'^ ^) is proportional to 
\5\. From the von Schweidler law we then arrive at 



AC CX |cr 



-l/(2a)-l/(2fc)|^|l/b+l 



\a\ < \5\ < 1. (42) 



Hence, the exponent of von Schweidler's law is also 
present in the behavior of AC, as can bee seen in Fig. 11 
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where A(C) is plotted as a function of 1/|5| for forces be- 
low the delocahzaton threshold, in the liquid. Note that 
for i^ox ^ 0, i.e., (5 — > — 1, the von Schweidler regime is 
cut off due to the microscopic relaxation of the correla- 
tion functions, and A({Fcx — 0) — ~ F^^ is obtained 
as mentioned further above. As the glass transition is ap- 
proached, this window of initial linear response shrinks. 



Again, Eq. (42) includes a scaling prediction: plotting 



AC as a function of S, where the rescaled variables are 
defined above, the curves for all e close to the glass tran- 
sition can be scaled onto one master curve. This is shown 
in the inset of Fig. [Tl] The von Schweidler law is seen for 
(5 <C 1, in the form A^ ~ (5^^''+-'^. This exponent is close 
to 2.6 for the model parameters chosen here, as one can 
verify in the figure. 



C. Large-Force Plateau 

We come briefly back to the case, where the probe par- 
ticle is delocalized. For S ^ 0{1), the /3-scaling regime 
does not give the dominant contribution to A^ any longer 
and the short-time dynamics of the correlation functions 
is more important. In this case it is sufRcient to consider 
the probe dynamics alone and that of the host liquid as 
essentially arrested, i.e., Vj(j){t) sa Vj^^^ as constants. The 
friction coefficient for S — > oo, ACoo — ^Ci^cx ^ oo) can 
then be calculated as 



ACoo^^oy {(l-M)Re0f(i)+/i01(t)} 

-0o{(l-M)Im(^^(z = O)+Mlm,^^(z = O)} . (43) 



From Eq. (6b), tp'^ — (j) [z — Q) fulfills a linear equation 

system, = i^^Jr^4>Q + (± ['4''^ : 4'q] ■ This yields an alge- 
braic approximation for ACoo- For the model discussed 
here, we get 



ACoo^M/cff + o(i/^;l). 



(44) 



where fcs — 0{1) is the approximately constant value of 
the host-liquid correlator (j){t) for the time window over 
which the probe-particle correlations decay. The incre- 
ment of the friction coefficient as i^ex oo hence is given 
by the parameter /x that controls the admixture of "per- 
pendicular" modes to the integral determining A^. This 
fact has been exploited in Ref. 10 in order to quanti- 
tatively describe the large-force amplitude of measured 
data. Setting fi — results in an expression for A^ that 
decays to zero as l/i^ox large forces. Setting further 
ti| = 0, one obtains the model used in Ref. ,9!, where 



ACoo = /cff(l + ^^e'ff)/(l 



^cff 



This expression 



holds for v^g > 1, ensuring positivity of the result. 

The reason that the correlation function (t) does not 
contribute to the friction-coefficient increment at large 
forces, is its oscillatory decay. Fig. [12] shows exemplary 
results for a very large force F^yr- both (j){t), shown as dot- 
ted lines, and (/'^^(i) decay roughly exponentially, with a 
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FIG. 12. Correlations functions <^jj(t) and (j>±{t) of the 
schematic model, with parameters as in Fig. [l] but a large 
external force, _Fcx = 402.7. Full lines show the numerical 
solutions of the model, dashed lines the least-square fits using 
an exponentially damped function <?!>j|(i) = e~''*e*'^°'="'"*'* with 
parameters v — 22.6 and ^ — —2.1. The host-liquid correlator 
<f){t) is shown as a dotted line in the top and bottom panels. 



relaxation time that is much larger than the one relevant 
for 0||(t) in this regime. The latter becomes proportional 
to a damped oscillation, (/'jj(t) « exp[— i + in^^Fcxt] in 
the absence of any memory- kernel damping, cf. Eq. ( |3a[ ). 
Taking into account the remaining damping from the 
memory kernel, we find that (/)| (t) for large i^ox is nicely 
fitted with exp[— z^i] exp[i(K||i\,x + OA- Such a fit is in- 
cluded in Fig. [12] with parameters i/ = 22.6 and ^ = —2.1 
determined by a least-square error minimization; it is 
virtually indistinguishable from the numerical result for 
(/)| (t) in the figure. 

The microscopic interpretation of these oscillations is 
a steady motion of the probe, r^{t) — r'*(0) +vt, leading 
to (j)g{t) = exp[iq ■ vt]. For q \\ Fg^, choosing a wave vec- 
tor characterizing typical probe-host interaction length 
scales, we recover ({Fcx-^ oo) = The schematic- 

model parameter Ky therefore relates typical length scales 
to the high- force friction coefficient. It is also clearly seen 
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that for q _L Fgx, no oscillations remain in the correla- 
tion function; this is correctly captured in the schematic 
model. 



VII. CONCLUSION 

We have discussed a schematic mode-coupling model 
for the nonlinear response of the friction coefficient in 
force-driven active microrheology. Expressions describ- 
ing the asymptotic behavior of the probe-particle den- 
sity correlation function in the directions parallel and 
perpendicular to the force within the schematic model 
have been derived. From these, together with the known 
asymptotic behavior of the equilibrium host-liquid corre- 
lation function, we have inferred two-parameter scaling 
laws that yield scaling forms for the probe-particle fric- 
tion coefficient. The two small parameters are the dis- 
tance to the glass transition of the host liquid, e, and 
the distance S to the delocalization transition, where a 
frozen-in probe becomes mobile in response to a force- 
induced local melting of the host. There exists a scaling 
limit based on the combination |£/(5^|: if this parameter 
is small, the external force is a small perturbation to the 
glassy dynamics of the host liquid, while for \S\ <^ ^/lej, 
the vicinity to the force-induced local melting of the glass 
dominates the dynamics. 

The asymptotic expressions for the correlation func- 
tions lead to the identification of five relevant time scales: 
two scales, to- and t'^, determine the power-law regimes 
for the host-liquid correlator. These diverge upon ap- 
proaching the glass transition, cr — > 0, and are indepen- 
dent of the external force characterized by 6. For the cou- 
pling of the probe correlators to that of the host liquid, 
a time scale ts was identified that diverges as the critical 
force is approached, (5—^0. It separates two regimes of 
proportionality, G^{t) cx G{t), but with different prefac- 
tors. Considering the long-time decay of the probe corre- 
lator, we identified two further time scales, ti/2 and t'^g, 
describing, respectively, the decay inside a glassy host 
liquid and the coupling of the probe relaxation to the 
host's von Schweidler law. For the friction coefficient, 
the various asymptotic regimes of the correlation func- 
tion induce several power-law regimes around the critical 
threshold force, F^^, marking the point where the exter- 
nal force becomes more effective in breaking cages than 
thermal fluctuations. 

In the linear response regime, the well-known asymp- 
totic result of equilibrium MCT is recovered, stating that 
the structural relaxation of the probe is, up to a wave- 
vector dependent amplitude, given by the relaxation of 
the host, cf. Eq. (20). In a sense, this is the justification 



typically of the size expected from the Stokes- Einstein re- 
lation |34j- Note, however, that the true Stokes-Einstein 
relation is of the form D ^ k^T/rj, where the factor 
k^T identifies it as a relation connecting a macroscopic 
transport coefficient to microscopic fluctuations. Indeed, 
recent experiments indicate that the factor k^T is absent 
in the empirical relation between D in rj found in viscous 
liquids [35^ , in agreement with Eq. ( 20 1 . 



Small external forces will preserve this coupling, but 
as the critical force is approached, the range of validity 
of Eq. (20) continuously shrinks. The initial deviation 



of the discussion of microrheology response in terms of 
generalized Stokes-Einstein relations, where one tries to 
relate the diffusivity of the probe particle to the collective 
dynamics of the surrounding host (expressed by its shear 
viscosity). Within MCT, a relation of the form D ^ l/ry 
is the consequence of Eq. ( 20 ) , with prefactors that are 



from the linear reponse regime, by symmetry, has to be 
quadratic, A({Fex) — A({Fcx = 0) ~ —F^^. Approaching 
F^^, however, the von Schweidler law governs the force 
thinning behavior, i.e., a power law involving the expo- 
nent b is seen in A^. This is, however, a small effect, not 
immediately apparent from experimental data. 

The delocalization transition is, within the present 
schematic model, identified as a continuous transition, 
i.e. the nonergodicity parameter of the probe-particle 
density correlators vanishes continuously. In a micro- 
scopic model, this would correspond to a continuous 
broadening of the probability distribution function for 
the probe position. This can be understood as a single- 
particle localization length that diverges due to the ap- 
plied force. However, in the schematic model, no such 
interpretation is obvious, since the model does not carry 
any information related to large length scales. It is also 
expected that the exponents found in the asymptotic 
analysis carried out here will be modified by taking into 
account long- wavelength fluctuations (50] . 

Above the threshold force, tracer motion is always de- 
localized. For this case, the friction coefficient increment 
AC follows two power laws: a trivial force thinning with 
exponent —1 is predicted in the glass, valid for very small 
I Fox — ^^cxl- This law is, closer to the glass transition, re- 
placed by a power law involving the MCT critical expo- 
nent, i.e., force thinning with exponent y = 1 — 1/a < — 1. 
Inserting typical values found for hard-sphere like sys- 
tems within MCT, a ss 1/3, one gets for the force thin- 
ning exponent y w — 2. 

Considering still larger forces, F^x — > oo, a feature seen 
in experiment is a second plateau in A^ > 0, indicating 
that the friction experienced by the pulled particle is not 
just the solvent friction. Although hydrodynamic inter- 
actions (HI) in the colloidal suspension will play a major 
role in this regime ^32j , the second plateau is not only due 
to them, as it is also found in Brownian-dynamics simu- 
lations that do not include HI. Within the schematic 
model, the additional contribution can be traced back 
to the correlation function 4'']_{t) mimicking the density 
correlators in the direction perpendicular to the force. 

Although the schematic model we have discussed gives 
a reasonable quantiative description of the friction co- 
efficients, and a qualitative one for intermediate-length- 
scale density correlation functions, the question remains 
whether the force-driven delocalization of a probe par- 
ticle in the glass is continuous in the sense mentioned 
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above. Another possibility is that of a discontinuous 
local yielding, where the plateau of the tagged-particle 
correlation function does not decrease to zero as 
Fex — >■ fex ~ Oj but to somc nonzcro constant (so that 
a discontinuous jump in the nonergodicity parameter at 
F^^ results). These scenarios are not easily distinguished 
following the behavior of the friction coefficient. Corre- 
lation functions from computer simulation of Brownian 
soft-sphere systems Unj appear compatible with the 
continuous scenario incorporated in the schematic model 
discussed here. However, recent MD simulations of the 
mean-squared displacement in active microrheology show 
a plateau that does not change appreciably with varying 
-Fex [S], and thus may point to a discontinuous transi- 
tion scenario. Note that even the connection between 
schematic and "full" microscopic MCT remains unclear 
at the moment. This connection hinges upon the assump- 
tion that the two modes chosen for the probe correlator 
in the schematic model are sufficient to represent the full 
set of critical modes in the bifurcation transition of the 



full model. In this respect, the microscopic model of 
Ref. is an extension of MCT that is qualitatively dif- 
ferent from other extensions of the original MCT where 
it could be shown that the latter's bifurcation class Ai 
was not left. 
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